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ABSTRACT 

This paper presents calculations for forbidden emission line profile shapes arising from 
colliding wind binaries. The main application is for systems involving a Wolf-Rayet (WR) 
star and an OB star companion. The WR wind is assumed to dominate the forbidden line 
emission. The colliding wind interaction is treated as an archimedean spiral with an inner 
boundary. Under the assumptions of the model, the major findings are as follows, (a) The 
redistribution of the WR wind as a result of the wind collision is not flux conservative but 
typically produces an excess of line emission; however, this excess is modest at around the 
10% level, (b) Deviations from a flat-top profile shape for a spherical wind are greatest for 
viewing inclinations that are more nearly face-on to the orbital plane. At intermediate viewing 
inclinations, profiles display only mild deviations from a fiat-top shape, (c) The profile shape 
can be used to constrain the colliding wind bow shock opening angle, (d) Structure in the line 
profile tends to be suppressed in binaries of shorter periods, (e) Obtaining data for multiple 
forbidden lines is important since different lines probe different characteristic radial scales. 
Our models are discussed in relation to Infrared Space Observatory data for WR 1 47 and 7 Vel 
(WR 1 1). The lines for WR 147 are probably not accurate enough to draw firm conclusions. 
For 7 Vel, individual line morphologies are broadly reproducible but not simultaneously so 
for the claimed wind and orbital parameters. Overall, the effort demonstrates how lines that 
are sensitive to the large-scale wind can help to deduce binary system properties and provide 
new tests of numerical simulations. 

Key words: stars: binaries - stars: early-type - stars: mass-loss - stars: winds, outflows - 
stars: Wolf-Rayet 



1 INTRODUCTION 

Massive star colliding wind systems continue to be an area of in- 
tense study. The colliding wind interaction region (CWIR) is intrin- 
sically interesting as a problem of hydrodynamic flow (Girard & 
Willson 1987; Shore & Brown 1988; Stevens, Blondin, & Pollock 
1992; Usov 1992; Canto, Raga, & Wilkin 1996; Gayley, Owocki, 
& Cranmer 1997). The CWIRs lead to a number of interesting ob- 
servable effects as well, such as the X-ray emissions (e.g., Stevens 
etal. 1996; Maeda et al. 1999; Zhekov & Skinner 2000; Pittard 
et al. 2002; Henley, Stevens, & Pittard 2003), ultraviolet (UV), op- 
tical, and infrared (IR) emission line profile effects (e.g., Stevens 
1993; Liihrs 1997; Stevens & Howarth 1999; Hill et al. 2000; Hill, 
Moffat, & St-Louis 2002), dust emission (e.g., Tuthill, Monnier, 
& Danchi 1999; Monnier, Tuthill, & Danchi 1999; Tuthill et al. 
2006; Tuthill et al. 2008), and polarization (e.g., Brown, McLean, 
& Emslie 1978; Drissen et al. 1986; Villar-Sbaffi et al. 2003, 2005, 
2006). Of chief interest between relating the observables to the hy- 
drodynamical predictions are the determinations of the mass-loss 
rates for the two stars, their orbital properties, and the opening an- 
gle of the colliding wind bow shock. Obtaining the orbital solution 



implies deriving the masses of the two stars, which is important 
for understanding massive star evolution. As is typical with unre- 
solved binary systems, the viewing inclination i creates a challenge 
to obtaining the full orbital solution and stellar masses. 

In the study of CWIRs of massive binaries, the interpretation 
of optically thin lines remains prominent, such as X-ray lines that 
encode information about opening angle and viewing inclination. 
One area that has been largely ignored is the influence of the CWIR 
on the shape of forbidden emission line profiles, especially those 
in the IR. The Infrared Space Observatory ('ISO') and the Spitzer 
Space Telescope ('Spitzer) have provided a large quantity of high 
quality IR spectra for massive stars. Especially in the case of the 
evolved Wolf-Rayet (WR) stars, the IR spectra are rich in emission 
lines and in particular forbidden lines. These forbidden lines have 
been useful for deriving or constraining gaseous abundances that 
are relevant for testing models of massive star evolution (e.g., Bar- 
low, Roche, & Aitken 1988; Willis et al. 1997; Dessart et al. 2000; 
Morris et al. 2000; Morris, Crowther, & Houck 2004 ; Smith & 
Houck 2005; Ignace et al. 2007). 

An interesting property of these forbidden lines is that they 
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form at large radius in the wind of a WR star. Consequently, the 
line formation occurs well beyond the wind acceleration thus sam- 
pling the flow at the wind terminal speed Voa. Being optically thin 
these forbidden lines should be flat-topped in shape, and so they are 
excellent for measuring Voc values (Barlow et al. 1988). 

But the flat-topped profile morphology only holds if the wind 
is spherically symmetric. If a WR star is part of a binary, then even 
if the WR wind is intrinsically spherically symmetric, the CWIR 
ensures that it will not remain so. Equally interesting is the fact 
that forbidden lines can form over tens of thousands of WR stellar 
radii (e.g., Ignace et al. 2001), corresponding to scales of several 
AU and comparable to binary separations in massive star binaries. 
As a result, forbidden emission profiles should not be flat-topped 
in general in these systems, and their deviation from flat-top - like 
the case for other CWIR emission line diagnostics - can be used to 
constrain the properties of the CWIR geometry. 

In this paper we calculate a range of forbidden line pro- 
file shapes as a function of wind and orbit parameters. Section 2 
presents a brief overview of forbidden line formation, a derivation 
of emission profiles in a spherical wind, and then in a colliding 
wind system. Then Section 3 describes applications to two systems: 
WR 147 and 7 Vel. Concluding remarks are given in Section 4. 



2 FORBIDDEN LINE FORMATION IN WINDS 

Our presentation of the theory of forbidden line emission gener- 
ally follows Osterbrock (1989) regarding the atomic physics, Bar- 
low et al. (1988) for application to winds, and Ignace & Brimeyer 
(2006) for notation. The following sections describe (a) the two- 
level approximation for forbidden lines, (b) the solution for the for- 
bidden line emission in a spherical wind, (c) the adopted model 
used to describe a CWIR between a WR star and an OB star, and (d) 
adjustments to the forbidden line profile calculation arising from 
the CWIR. 



2.1 The Atomic Physics 

For simplicity the two-level atom approximation is adopted for a 
fine structure transition of an ion species. The lower level will be 
' 1 ' and the upper level '2' . It is common to introduce a critical den- 
sity n c that signifies the transition from the higher density medium 
where de-excitations are dominated by collisions versus the lower 
density zone where spontaneous decay dominates. Since excitation 
is assumed to derive from collisions only, the emissivity transitions 
from a function that is linear in density in the collisional regime 
to one that is quadratic in density in the decay regime. The critical 
density is defined as 

A21 



where 



n c = 



q>i 



(1) 



where A21 [s ] is the Einstein A-value for the transition, and 521 
is the downward volume collisional de-excitation rate [cm 3 s -1 ]. 
The volume emissivity j [erg s _1 cm -3 sr" 1 ] is 



1 L 

J = — I1V21 A21 m- 

4-7T 



(2) 



In the two level atom, the entire density of elemental species E in 
ion stage i is given by tie,i = ni + "-2- One can derive the ratio 
ri2/ni from equilibrium conditions and use that to solve for the 
population of the upper level: (Barlow et al. 1988): 

th,b n e qi2 



n 2 



n e qi2 + n c q 2 i + A21 



(3) 



*? 21 _ 91 
912 32 



with 



(3 = hv 2 i/kT e 



(4) 



(5) 



for V21 the frequency of the line transition, T a the electron temper- 
ature, and gj the statistical weight of the level, with gj — 2 J + 1. 
The downward collisional volume rate is given by 



q>i 



8.629 x 10" fa U12 



1/2 



92 



(6) 



where W12 is the collision strength. A summary of line transition 
data and critical densities is given in Table □ 

Forbidden emission lines in stellar winds are optically thin. 
The above relations are typically combined with the intent of inte- 
grating the emissivity over the wind and solving for the ionic abun- 
dance, derivable from observed forbidden emission lines and used 
to constrain gas abundances in evolved stars to test massive star 
evolution models. The primary goal of our paper is quite different. 
We seek to derive line profile shapes to determine to what extent 
those shapes may be used to deduce the properties of binary or- 
bits and stellar winds in colliding wind systems. To do so, we first 
review the solution for a spherical wind. 



2.2 Spherical Winds 

It is useful to review briefly the line emission for a forbidden line 
from a spherically symmetric wind, both to establish notation and 
as a reference to use with the non-spherical case. For a spherical 
wind with mass-loss rate M and radial speed v(r) for radius r, the 
mass density is 



Psph 



M 



4-7T r 2 v(r) 



(7) 



The electron number density is n c = p sp h/ fJ. e mH, for fi c the mean 
molecular weight per free electron. 

The region of line formation is set roughly by the radius at 
which rt c = n c . For hot stars M ~ 10~ M© yr _1 and larger, 
with Voo ~ 1000 km s _1 or more for OB stars (e.g., Lamers & 
Cassinelli 1999). To estimate the radius of line formation, we in- 
troduce an electron number density scale factor that depends on 
these basic wind parameters: 



no,e = 



M 



4ir Rl Voc fJ-c run 
Then the critical radius becomes 



R* 



np,» 



(8) 



(9) 



With no,e of order 10 13 cm 3 for a WR star, and n c about 10 5 
cm" 3 , forbidden line emission forms in the far wind, around 
10 4 i?». 

To determine the emission line shape, we combine equa- 
tions l[2} and l[3} to obtain the emissivity: 



J — Jo 



D c n e /no,, 



1 + uj + uj D c n e /no,, 



(10) 



where jo is a constant that depends on the line of interest, D c ^ 
1 is the 'clumping factor' that introduces a dependence on the 
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Table 1. Forbidden Line Data" 



Ion 


Transition 


9i 92 A 


P 


VJ 


A 2 i 


n b 








Om) 






(10" 3 ) s" 1 




(10 4 )cm" 3 



CalV 


2p , 
*l/2 " 


~* ^3/2 


4 


2 


3.2 


0.45 


3.14 


545 


1.06 


1200 


Nell 


2 Pl/2 - 


- 2 P 3 /2 


4 


2 


12.8 


0.113 


2.24 


8.55 


0.30 


65.0 


Nelll 


3 Pl- 


- 3 P 2 


5 


3 


15.6 


0.092 


1.83 


5.97 


1.65 


12.6 


SlV 


2 P 3 /2 " 


^ 2 Pl/2 


2 


4 


10.5 


0.137 


0.57 


7.73 


6.42 


5.6 


Oiv 


4 P 5 /2 - 


- 4 Pl/2 


2 


6 


32.6 


0.044 


0.35 


0.518 


0.69 


5.2 


NeV 


3 P2- 


-. 3 Pl 


3 


5 


14.3 


0.101 


0.66 


4.59 


5.82 


4.6 


Sill 


3 P 2 - 


-. 3 Pl 


3 


5 


18.7 


0.077 


0.65 


2.07 


5.81 


2.1 


Sim 


3 P 2 - 


^ 3 Pi 


3 


5 


38.2 


0.038 


0.62 


0.242 


10.4 


0.13 



a Atomic data taken from Pradhan & Peng 1995. b Evaluated for T c = 10 4 K. 



wind clumping (following Dessart et al. 2000). Note that Ignace & 
Brimeyer (2006) did not account for clumping in their expressions. 

The total luminosity of the optically thin line emission is given 
by a volume integration involving the emissivity: 



Li — 4-7T / j(r) D c 1 r 2 dr d cos § dip. 



(11) 



where (r, 1?, ip) are spherical coordinates in the star system, and the 
inverse of the clumping factor is the volume filling factor. Substitut- 
ing for a normalized radius x = r/R t and inserting the emissivity 
function above, one has 



dx d cos 1? dtp 



Lo 

4tt J 1 + vj + wD~ x (x/xc) 2 ' 



(12) 



where Lo = 167r 2 jo.R 3 , and the line is assumed to form predomi- 
nantly over the constant velocity flow where par -2 since x c 1 
(c.f., Fig. 1 of Ignace & Brimeyer 2006). 

The angular quantities integrate to 4-7T, and the radial integral 
is of a standard form. Given that x c 2> 1, the total line luminosity 
becomes 



L - 



n/2 



(13) 



Note that a constant clumping factor appears as a square root coef- 
ficient to the critical radius, and so the product acts essentially as 
a 'transformed' critical radius. Also note that the scale constant Lq 
depends on the mass-loss rate. The assumption is that equation d 1 3 b 
holds if M is already a dumping-corrected value. 

This of course is the total volume-integrated emission. More 
relevant to our study is the emission profile shape. Optically thin 
lines that form in a constant expansion and spherically symmetric 
wind produce flat-topped emission profiles. The isovelocity zones 
are given by 



"Woo Mi 



(14) 



where 8 is the polar angle from the observer axis z. With v z fixed, 
the isovelocity zones are cones. Introducing w z — v z /voo, and 
noting that dw z — —d/j,, the profile shape dLi/dv z = A -1 L v 
becomes 



dLi 
dv z 



lLi_ 

2 Vao 



(15) 



In the consideration of line shapes from colliding wind systems, the 
above constructions will continue to prove useful. 



WR 



. 2tt — 9? 



r(<p) 



CWIR 



Figure 1. Schematic of the spiral Colliding Wind Interaction Region 
(CWIR) in the binary orbital plane. The WR star is at left and the OB com- 
panion at right. In this "snapshot" the binary separation is r or b . The orbital 
motion is counterclockwise as signified by the angular speed lo, and so the 
spiral CWIR winds out in a clockwise fashion. A point along the CWIR P 
is indicated as having polar coordinates r and 2tt — ip. 



2.3 Colliding Wind Binaries 

Ignace & Brimeyer (2006) explored the effects of axisymmetric 
density distributions for forbidden lines shapes from single star 
winds. In that study the key parameters that determined deviations 
from flat-top morphologies were the viewing inclination, the lati- 
tudinal density description (whether bipolar or disk-like flow), and 
the ionization balance. The colliding wind systems allow for con- 
siderably more influences, in part because the system is completely 
non-axisymmetric and also because there are really two sets of pa- 
rameters to consider - those of the orbit and those of the stellar 
winds. 

In order to explore the range of line shapes and their relevance 
for observations, some approximations will be adopted. We con- 
sider WR+OB systems such that it is the WR wind that dominates 
the emission of the forbidden line. This is reasonable because the 
mass-loss is much higher and the abundances are non-solar for the 
WR component. For example, [Nelll] at 15.56 microns is typically 
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prominent across WR subtypes (Ignace et al. 2007), and neon can 
be enhanced above solar by factors of a few or much more in a 
WC wind. As a result, one major assumption of our study is that 
the confined wind of the OB companion as dominated by the more 
massive WR wind will have zero emissivity. 

As noted before, the forbidden lines form at large radii (i.e., 
compared to the stellar sizes, but not necessarily larger than the or- 
bit size). We will assume that the WR wind is in constant expansion 
for evaluating the line profile shape. 

For the binary we choose to work in a frame centered on the 
WR star (see Fig. [B. The OB companion is taken to follow an el- 
liptical orbit of semi-major axis a and eccentricity e. The observer 
is located on the z axis, and the normal to the orbit plane is defined 
to be z, . The viewing inclination is given by cos i — z ■ z*. 

The WR coordinate system is described by spherical coor- 
dinates (r,i9,(/p), and Cartesian coordinates (x*,y*, z*). The ob- 
server frame also centered on the WR star has coordinates (r, 9, a) 
and (x,y,z). The x-z and x*-z* planes are taken as coincident. 
Azimuthal angles are measured counterclockwise in standard fash- 
ion from the x and x* axes. The reference point for the orbit is 
the periastron location ip p . At any given epoch, the OB star will be 
located at azimuth tp s , which we take to be the orbital phase. 

The CWIR will be described in a largely asymptotic form. 
We consider the interaction region to have two primary zones: a 
shocked layer for the OB star wind, and a shocked layer for the 
WR wind. The interface is taken to describe the confinement of the 
OB wind. Since we ignore any line emission from the OB wind, 
the trace of that interface and its interior is treated as an emission 
'cavity' . Line emission then arises from two sectors only. The first 
is the shocked WR wind. The second is an approximately spherical 
WR wind that exists outside the shock layer. 

The asymptotic shape of the CWIR is that of a spiral, mo- 
tivated by observations of 'pinwheel' nebula from dust and radio 
emission in several WR binaries (e.g., Monnier et al. 2002). If the 
OB companion is at azimuth ip a and orbiting counterclockwise, 
then the interaction region spirals outward in a clockwise direction, 
and vice versa. We assume that the bowshock at the OB star has 
circular cross-section. The 'memory' of that interaction is moving 
radially outward at constant speed. The intersection of the spiral in- 
teraction region with a spherical shell centered on the WR star is as- 
sumed to maintain a constant solid angle of circular cross-section. 

The following relations describe the CWIR as related to the 
binary orbit. The OB star is located at a radius r or b in the orbital 
plane as given by 



^"orb 



1 + e cos(<p — ip p ) 



(16) 



At this distance the interface between the WR and OB star winds is 
taken to have a half-opening angle (3. The shocked layer of the WR 
wind is taken to have a half-opening angle f3' . Thus the intersection 
of a spherical shell with the compressed WR wind layer appears 
as an annulus of angular width A/3 = /?' — (3. The characteristic 
lateral lengths I and I' associated with these angles are: 



r orb tan/3, 
r orb tan/3'. 



(17) 
(18) 



Assuming the shock is strong, the immediate post-shock density 
will be 4 times that in the normal spherical WR wind component. 
Assuming an adiabatic shock layer with constant post-shock den- 
sity, then the annular cross-section of the shocked layer has one 



quarter the solid angle of the confined OB wind, leading to the fol- 
low relation that we adopt for our models: 



cos 



5 cos (3 ■ 



(19) 



Note that we treat f3 largely as a free parameter in our models 
to explore profile shape effects. In fact, its value should be related 
to the ratio of wind momenta for the WR and OB star winds. We 
define 



JOB 



(20) 



Then the stagnation point between the two stars will be located at 
a radius ro from the WR star as given by (e.g., Shore & Brown 
1988): 



ro 



^*orb 



(21) 



Recently, Gayley, Parsons, & Owocki (in preparation) have consid- 
ered the opening angle associated with a purely adiabatic shock, 
which is the case adopted for our work, in contrast to a radiative 
shock (see Canto et al. 1996; Antokhin, Owocki, & Brown 2004). 
Gayley et al. (in preparation) derive the following transcendental 
relation for the shock cone opening angle: 



cos /3 = 



2c 



1 



(22) 



This expression will be used in our application to WR 147 later in 
this paper. 

Interior to the orbit, the WR wind is assumed spherical. The 
bow shock region in the vicinity of the OB star will be dominated 
by hot, X-ray emitting gas (e.g., Parkin & Pittard 2008). Conse- 
quently, the requisite ions for contributing to the forbidden lines 
simply won't exist there. As a result, we model the bow shock 
zone OB star with cut-offs. We assume that the confined OB wind 
(hereafter referred to as the 'cavity') and the compressed WR wind 
(hereafter simply the 'compressed layer') extends to the stagnation 
point ro between the two stas and is too hot to produce forbid- 
den line emission for the ions of interest. We further assume that 
the gas remains hot some distance downstream of r or b. We adopt 
8ro for this distance, using Cassinelli et al. (2008) for adiabatic 
shocked wind flow around a blunt object as a crude guide. As long 
as a < r c or o > r c , the details of the bowshock region will 
not greatly impact the forbidden line profile shape. It is more rele- 
vant when a ~ r c , in which case our results can only be taken as 
illustrative. 

It remains then only to trace the center of the cavity with dis- 
tance and azimuth in the orbital plane to completely define the spi- 
ral interaction region through the WR wind. A cross-section of the 
spiral pattern is basically a 'shadow' of the wind collision that ad- 
vances through the WR wind at constant radial expansion. We sim- 
ply need an expression for the star's location around the orbit with 
phase to determine an equation of motion for this 'shadow' . Such a 
relation is derivable from Kepler's laws. 

At any given phase, the cross-section center advances radially 
according to 

r = v 00 . (23) 

Conservation of angular momentum C provides a relation for the 
star location: 



C 

1 '■ 

r orb 



(24) 
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Figure 2. Geometry related to the conical interaction region and the evalu- 
ation of eq. (28). The z point is toward the oberver and the z„ is the binary 
orbital axis. Point P signifies the axis of the conical CWIR, and point P' 
is located on the shock contact discontinuity. A similar set of spherical tri- 
angles would apply for a point on the boundary of the compressed layer at 
angles /3' and 9' . 

with the specific angular momentum given by 

C = y/GMa{l-e 2 ), (25) 

for M the summed mass of the two stars. Defining an angular ve- 
locity uj — 2n/P for period P, using Kepler's third law, and com- 
bining the three preceding equations, a differential equation for the 
cross-section center can be derived: 

^i = _^Zrkyi^, (26, 

ar r w 

where a convenient scaling parameter we call the 'wrapping length' 
r w is introduced: 

r w = Voa/oJ- (27) 

The wrapping length is related to the pitch angle of the spiral shape. 
Since r c is a scale for the emission of the line, it is natural to charac- 
terize models by the ratio 2nr w /r c — P/t c , where t c — r c /voo in- 
dicating the number of wrappings per critical radius crossing time. 
If the wind flow time across the critical radius is much longer than 
P, the spiral will have many circuits over that scale; but if the or- 
bital period is long, the spiral is essentially a cone over the span of 
a critical radius. 

2.4 Conical Interaction Regions 

Before tackling the case of general orbits, it is instructive to con- 
sider the extremely long period binaries for which the bow shock 
geometry is asymptotically conical in shape, having no spiral cur- 
vature. We will refer to this limit as the 'conical' (or 'linear') bow 
shock. 

In the approximation of negligible wrapping of the spiral, the 
cavity and compressed layer are cones centered on the line of cen- 
ters for the two stars. The isovelocity zones are also cones. An im- 
portant point is that we assume the flow in the compressed layer is 



Figure 3. Model line profiles in the 'conical bow shock' (i.e., r w S> r c ). 
Each profile is normalized to have unit area, and the horizontal axis is nor- 
malized velocity shift. Table[2]details the model parameters for each panel. 



radial and has the same speed as the WR terminal speed. This is a 
reasonable approximation for an adiabatic shock and a strong WR 
wind (see Tuthill et al. 2008). These isovelocity cones are centered 
on the observer's line-of-sight, and so are inclined to that of the 
CWIR. What is interesting is that the intersection of the interaction 
region with the isovelocity cones is a fixed pattern with radius. 

Consider a spherical shell. The cross-section of the isovelocity 
zones are rings. That of the interaction region is a ring also which, 
without loss of generality, is assumed to lie in the x-z plane. It is 
straightforward to find the crossing points between the two rings in 
terms of the azimuthal angle oto for the observer (see Fig. [2j . If the 
cavity center is at 9 C — 90° — i, the implicit solution for ao is 

cos (3 = cos 9 C cos 9 + sin f? c sin 6 cos ao (28) 

= w z w c + cos a v(l - - tof), (29) 

w z is the observed normalized velocity shift in the line and w c is 
the value for the isovelocity zone for the cavity axis. The same 
relation can be used for /3' and a' to determine the crossing points 
between an isovelocity ring on the shell and the circular boundary 
of the compression layer. The key is that ao and a' are constants 
with radius. 

There are several special cases that are notable. If viewed pole- 
on, intersections only occur for 9 > 90° — and 9 < 90° + /3'. 
When viewed edge-on, the axis of the interaction region coincides 
with the z-axis, in which case there is no emission for 9 < (3, and 
emission from the compressed layer is exclusively from f3 < 9 < 
P'. 

For emission from the compression region, the density is en- 
hanced by a constant factor of 4 in our treatment. In terms of the 
emission per unit solid angle, dL^/dQ., the integration along a 
hypothetical radial that lies entirely within the compression layer 
will be 8 times greater than for one in a purely spherical wind. Of 
course, the compression layer has a solid angle extent that is 1/4 
that of the cavity. In effect, the wind collision redistributes mass in 
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Table 2. Model Runs for Figure[3] 



Panel 


a/r c 


p 


i 






(°) 


(°) 


(a) 


0.1 


60 





(b) 


0.1 


45 





(c) 


0.1 


30 





(d) 


0.1 


60 


90 


(e) 


0.1 


45 


90 


m 


0.1 


30 


90 


(g> 


10 


40 


30 


(h) 


1 


40 


30 


(i) 


0.01 


40 


30 



a sector of the WR flow, and there is a net gain in line emission 
above what that sector would have produced if it were undisturbed. 
The wind collision leads in essence to a globally stuctured clump. 
The important point is that one should not expect CWIRs to con- 
serve flux in the forbidden line as compared to a spherical wind. 

Also, the CWIR does not extend down to the photosphere of 
the WR star, but is significantly offset in radius, as determined by 
the location of the stagnation point on the line of centers between 
the two stars and the width of the compressed WR wind. This 
means there is a minimum radius to the CWIR, interior to which 
a spherical WR wind makes a simple flat-top contribution to the 
line profile. 

Example profiles for the conical bow shock approximation are 
shown in Figure [3] for different binary separations relative to the 
critical radius a/r c , cavity opening angles P, and viewing inclina- 
tions i. Model parameters for the different panels are provided in 
Table[2] The principle conclusions are that: (a) only a pole-on view 
to the orbit produces a symmetric profile, with a double-horned 
appearance, (b) an edge-on view produces one that is maximally 
lopsided, and (c) generally an asymmetric double-horned profile 
shape results whose appearance relates to the viewing perspective 
and orbital parameters. Note that these profiles have been gaussian 
smoothed to simulate limited spectral resolution. Given that typical 
WR winds have ~ 1000 — 3000 km s _1 , smoothing with a 
gaussian of HWHM Sv/voo =0.1 was adopted to match roughly 
the resolution of ISO's SWS06 instrument (de Graauw et al. 1996). 

Finally, all of the examples in Figure[3]have excess emission 
relative to a purely spherical wind by factors of 10-20%. It happens 
that the total line flux as a function of the opening angle and binary 
separation is derivable analytically. There are four basic zones. As 
previously noted for radii r < r or b, the WR wind is spherical and 
contributes a flat-top contribution to the profile. In the cavity sector, 
there is no contribution. Then there is the spherical zone and the 
compressed layer for r ^ r or b. Accounting for these zones, the 
solution for the line luminosity, relative to a pure spherical wind, is 
denoted by A and given by 



A = 1 + - (1 - cos 0) 



tan 



tan 



'fir) 



where 



-Dei no 



(30) 



(31) 



1 + w 

In the limit that /3 = 0, equation d30t reduces to unity, as it must 




-2-10 1 2 

log ( a r c / r ) 

Figure 4. A plot of the integrated excess line flux for a colliding wind 
system as compared to a single spherical wind. This is plotted against the 
ratio of the critical radius associated for the line relative to the separation of 
the WR star and OB companion. An additional dimensionless multiplica- 
tive factor, a, accounts for the wind clumping and atomic constants for the 
transition (see eq. 1311 ). 



because there is no CWIR. For two identical winds, the opening 
angle is P = 90° , and the maximum line flux for a given value of 
auo/u c — aT c /ro becomes: 

"la uq^ 



A = 



1 — — tan 1 ( -\ H tan 1 

4-7T V u c J IT 



(32) 



If the binary orbit is exceedingly large, A = 1 is again recovered 
because the CWIR is displaced to a location of irrelevance with 
r rb 3> r c . On the other hand, if r c 3> r or b, all the arctangent 
factors reduce ton/2, thus 

11 — 3 cos P 



A = 



(33) 



Formally this has a maximum of 1.75. Figure [4] shows (A — 1) as 
a percentage excess of line emission above the spherical value as a 
function of ar c /ro- The different curves are for different /3 values, 
from to 90 degrees in 10 degree increments, with the excess being 
greater for larger /3. The location of maximum percentage excess 
can be derived from equation i30i and occurs at ar c /ro = \/TT/2 
for all p. 

These curves are accurate only within the assumptions for the 
adopted geometry. A different set of curves would result for exam- 
ple if the relation between p' and P were different. Also, for ease 
of calculation, the bow shock exists for r > r ro b and not at all for 
interior radii, which is incorrect because this neglects the shape of 
the head of the bow shock. The exercise does serve to demonstrate 
that (a) colliding wind effects will generally not conserve line flux, 
(b) line flux excesses at the level of ~ 10% can be expected, (c) 
the excess will depend on the clumping factor in the large-scale 
wind, and (d) elemental abundance determinations for WR winds 
from forbidden lines formed in colliding wind systems are not dra- 
matically biased by the CWIR. However, it may be that changes in 
ionization between the nominal spherical wind and the compressed 
layer could have a larger systematic influence that would need to be 
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Figure 5. Shown are a series of model calculations to illustrate the range of profile shapes that result as a function of bow shock parameters, viewing 
perspective, and orbital parameters. The first three columns are for orbital periods of 1, 100, and 10,000 years, respectively. The fourth column is for a strictly 
conical bow shock (essentially an infinite orbital period, as in Fig.fSJ. From top to bottom, the viewing inclination is i = 0, 2, 10, 30, 45, 60, 80, 88, and 90 
degrees. The orbits are circular (e = 0). Notably, the profile shape is sensitive to the spiral structure, a reflection of the orbital motion, even for rather long 
orbital periods of a century; the effectively conical bow shock limit is not achieved until a period of about 10 millenia. 



investigated more carefully in detailed simulations. We attempt to 
include such affects in our consideration of applications to WR 147 
and 7 Vel in Section 3. 



2.5 Spiral Interaction Regions 

For spiral-shaped interaction regions when binary orbital effects are 
important, issues that influence the profile shape are similar as for 
the conical shaped regions of the last section. However, there are a 
few more important factors to bear in mind. The spiral is still sym- 
metric about the orbital plane. That means the line profile should be 
symmetric about zero velocity for a pole-on inclination of the sys- 
tem. However, when edge-on the spiral shape now allows the inter- 
action region to transit (eventually) every isovelocity zone, thereby 
tending to dilute the influence of non-sphericity in the profile shape. 
Another important factor is that the direction of orbital revolution 
with respect to the observer will be relevant. 
One expects the following trends to hold: 

(i) Larger opening angles lead to profile structure over a greater 
range of observed velocity shifts. 



(ii) Flat-top segments seen toward the line wings places a limit 

on i — f3' . 

(iii) The total line luminosity derivation of equation J30t for 
the conical interaction regions remains relevant for the spiral case. 
With the assumption of constant solid angles, the proportions of the 
cavity and compressed layer to a spherical shell are constant with 
radius. This means that a spiral interaction region is simply a coni- 
cal one with radius-depdendent phase lags. However, those curves 
require a correction for the extension of the emission cavity down 
to the stagnation point in regimes where that portion impacts the 
line emission. 

(iv) When the wrapping radius is relatively small, the profile 
shape tends to deviate less from a flat-topped appearance, for most 
viewing inclinations. 

(v) The compressed layer leads to line emission being sensitive 
to the outer spiral shape even when the wrapping radius is some- 
what large. 

The last point is interesting, and is highlighted in Figure [5] 
which also provides a test of our code for general orbits against 
the semi-analytic case of a conical bow shock. The figure is or- 
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Figure 6. A selection of profile shapes as various individual orbital parameters are varied. Labels indicate control parameters. Each profile is is plotted against 
normalized velocity shift, as in previous figures. Profiles are vertically offset for better display of profile shape effects. Specific parameter values are described 
in the text. 



ganized in four columns. From bottom to top, the viewing in- 
clination ranges from pole-on to edge-on through angles of i = 
0, 2, 10, 30, 45, 60, 80, 88, and 90°. From left to right for the first 
three columns, the orbital period is 1 year, 100 years, and 10,000 
years. Profiles in the far right column are for the purely conical 
CWIR described in the previous section. The orbits are taken as cir- 
cular, with r or b = 0.001r c so that the CWIR exists essentially ev- 
erywhere. The profile shapes have not been gaussian smoothed. The 
bow shock has /3 = 40° and the compressed layer has f3' = 45° 
in each instance. Note that the vertical scales (not shown) are the 
same within a row but vary between rows. 

For the longest period case, the profiles in the third column 
match those of the strict conical bow shock case in the fourth col- 
umn. What is surprising is that models with P — 100 years for 
which the wrapping radius is larger than the critical radius by two 
orders of magnitude deviate from the expectations of the conical 
bow shock case. The reason is that the emission in the compressed 
zone is very large. The factor of 4 increase in density combined 
with the density square emissivity dependence in the outer wind 
along with the narrowness of the zone leads to relatively 'spikey' 
features in velocity that appear in the line profile. Of course with in- 
strumental smearing, the significance of these features will be sub- 



stantially reduced. Also note that since the geometry is fixed, all 
the profiles should have the same flux, which we have confirmed. 

It is clear that diverse profiles can result. Even for the highly 
controlled examples of Figure [5] and with some instrumental 
smoothing, profiles can be largely symmetric or strikingly not so; 
they may have two peaks, three peaks, or possibly four; and some 
profiles can be nearly flat-topped in shape, or have flat-top seg- 
ments. 

It is difficult to explore the parameter space exhaustively ow- 
ing to the large number of free parameters, namely i, j3, r c /a, 
r c /r w , ip p , and ip, plus the orbit direction, so a limited selection 
of illustrative results are shown in Figure [6] These have been con- 
volved with a gaussian of HWHM Sv/voo — 0.1. Each panel is 
labeled. The two longer panels at right are for pole-on and edge- 
on viewing inclinations in which the semi-major axis takes values 
of a/r c = 0.01, 0.1, 1.0, five profiles for each case from bottom 
to top. The orbits are circular, and there are five opening angles 
P = 30, 40, 50, 60, 70 degrees. The orbital period is 1 year. 7 Vel. 

The four panels at left are for more restricted variations in the 
parameters as listed. Generally P — 1 year, a/r c — 0.1, e = 0.1, 
i — 60°, j3 = 40°, and the periastron and companion star are both 
located at ip p = ip s — for these models. The top panel shows 
variations in P from 0.01 to 100 years in factors of ten from bottom 



© 0000 RAS, MNRAS 000, 000-000 



Forbidden Lines from Colliding Winds 9 



Table 3. Wind and Binary Parameters for WR 147 



WR 147" 


WN8+BV 


M 


25 X 1CT 6 M yr- 1 


Voa 


950 km s" 1 


R» 


21 R Q 


Me 


3.1 


Del 


10 




1.2 X 10 11 cm" 3 


e 


1/90 


Orbit 6 


d 


630 pc 


P 


1000 years 


a 


870 AU 


e 





i 


65° 


Interaction Region 




40° 


Tw 


330,000 i?»/rad 


r w /r c (Caiv) 


3300 


r w /Vc (Nem) 


330 


r w /r c (SlV) 


220 



Morris et al. 2000. 6 Churchwell et al. 1992; orbital parameters are 
highly uncertain, but given the distance and apparent separation of the two 
stars, a ~ 370/ cos i AU. c Using Gayley et al. 2009. 

top. For eccentricity, e = 0, 0.05, 0.1, 0.3, 0.8, again from bottom 
to top. The periastron location varies from to 160 degrees in 40 
degree increments; the same variation is used for the location of the 
star. Note that relative to a spherical wind, the line fluxes are never 
in excess by more than 23% for this set of models, with an average 
excess of about 7%. 



3 DISCUSSION 

There are currently relatively few WR colliding wind systems that 
have adequate data quality for the application of our model to in- 
terpret the CWIR geometries. ISO has a spectral resolving power 
of about 3000 but observed relatively few WR stars, and except for 
the fastest WR winds, the spectral resolving power of Spitzer is too 
low to resolve the features that our model predicts. However, ISO 
did observe some colliding wind binaries in its highest resolution 
mode, and we use two of them as test case studies: WR 147 and 
7 Vel (WR 1 1). The first is a wide binary with an unknown period, 
although it is likely on the order of millenia, and the second is a 
relatively close binary with an orbital period of a few months. 

3.1 Application to WR 147 

Adopted parameters for the WR 147 system are provided in Ta- 
ble [3] The separation of stellar components is at least 360 AU, 
or about 10 5 cm, based on their projected separation in the sky 
(Churchwell et al. 1992). The period is estimated to be at least 
P ~ 10 3 years. The wrapping radius is thus about > 10 17 cm. The 
wind density scale is around no, c ~ 10 11 cm" 3 . We assume a sta- 
tistically likely viewing inclination of 65°, in which case the orbital 



separation becomes 870 AU. Having no other knowledge about the 
orbit, we adopt e = 0. Radio observations by Williams et al. (1997) 
shows that the bow shock head is quite close to the companion star, 
for which they argue for a BV companion. The ratio of wind mo- 
menta is then e ~ 1/90. Using equation J22t for the bow shock 
opening angle for an adiabatic shock, we obtain (3 « 40° . 

Lines of [Caiv], [Nem], and [SlV] were observed with ISO 
(Fig. 8 of Morris et al. 2000). The WR wind density at the distance 
of the B star is n c ~ 10 3 cm" 3 , well beyond the critical radii 
values for these three lines. Deviations in the emission lines from 
flat-top should be mild, consistent with the observations of Mor- 
ris et al.With the sulfur line forming farthest out, we might expect 
that some structure would be observed in it. In addition, the stellar 
components are separated by about 10 4 i?», whereas the winding 
radius is around 30 (or more) times that value, thus the CWIR is 
closely approximated by the conical bow shock model described in 
section [2~4l 

Synthetic spectral profiles are displayed in Figure [7] Each 
panel shows model profiles plotted against normalized Doppler 
shift. The different profiles are for different orbital phase location, 
sampling a full orbit in 30° intervals. The model profiles have been 
convolved with a gaussian of Av/voo = 0.2 to simulate ISO's 
spectral resolution for the relatively lower Voo wind of WR 147. 
The line of [Caiv] has the smallest critical radius and shows an 
essentially flat-top shape. The neon and sulfur lines show greater 
variation, as their critical radii are somewhat closer to the B star 
companion. The observations of Morris et al. (2000) are sugges- 
tive of a double-horned feature in [Nem] and [SlV], and indeed our 
models manage to produce such a profile shape in narrow ranges 
of orbital phase. However, it is unclear to what extent the profile 
shapes may be trusted. A re-analysis of the [Caiv] line by Ignace 
et al. (2001) suggests an asymmetric line shape. If real, that could 
indicate an elliptical orbit. A series of model runs with large eccen- 
tricites greater than 0.9 retain double-peaked morphologies in the 
Ne and S lines and shows some small deviations from flat-top in the 
Ca line for a range of orbital phases. Although better S/N data at 
high spectral resolution is needed to measure the line shapes more 
accurately, this preliminary application shows intriguing potential. 

3.2 Application to 7 Velorum (WR 11) 

Here we apply our synthetic line profile approach to the case of 
7 Vel consisting of a WC8 star with an 07.5 companion. Adopted 
parameters for the WR wind and the orbit are given in Table [4] 
This colliding wind system has been studied extensively. The star 
and wind parameters come from de Marco et al. (2000), and the or- 
bital parameters come from the interferometric study by North et al. 
(2007). Note that an independent study by Millour et al. (2007) 
yields mostly similar parameters. Four forbidden lines were ob- 
served with ISO: [Caiv], [Nell], [Nem], and [SlV] (see Dessart 
et al. 2000 and Ignace et al. 2001). In contrast to WR 147, the or- 
bital parameters are quite well-known, and so 7 Vel constitutes a 
significant test case for our modeling. We note that the adopted 
opening angle of f3 = 85° is significantly larger than the value for 
an adiabatic shock at « 45° (see eq. 11221 ) or a radiative shock at 
« 35° (see Canto et al. 1996). We use this larger value based on 
the X-ray study by Henley et al. (2005). 

There are a number of new issues that arise in the case of 7 Vel 
that are different from the wide binary WR 147, notably that (a) 
the winding radius is comparable to the critical radii of the ob- 
served forbidden lines and (b) the critical radii of most lines are 
only somewhat larger than the orbital semi-major (within a factor 
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Figure 7. Model line profile shapes as indicated with application to WR 147. Each profile in a column is for a different orbital phase location of the B star 
companion. 



of about 10), and for [Caiv], the critical radius is nearly equal to the 
separation of the stars at apastron. The three primary scale lengths 
of the problem - r w , r or b, and r c - are somewhat comparable for 
all four lines of interest. As a result, the single greatest deficiency 
of our model becomes potentially important, namely the treatment 
of the bow shock head. 

Our treatment of the shape of the bow shock interior to the 
instantaneous position of the companion star is extremely rough. 
We determine the stagnation point ro based on equation l !2 1 1 . The 
emissivity 'cavity' extends inward in a conical fashion from the 
companion star to the stagnation point. This is a gross misrepresen- 
tation of the rounded shape that the bow shock should take. How- 
ever, this sector of the bow shock is the portion where the two winds 
intersect most nearly head-on and which produces the observed X- 
ray emissions (e.g., Skinner et al. 2001; Henley et al. 2005). This 
means that the low ionization states of calcium, neon, and sulfur 
that contribute to the IR forbidden line emissions simply will not 
exist in this sector at all. 

Figure [8] left side, shows synthetic emission profiles for all 
four observed lines in a format similar to Figure [7] for WR 147. 
The model lines were gaussian smoothed with Av/voo = 0.14 
to simulate the resolution of ISO for the terminal speed of 7 Vel's 
wind. The profiles are only moderately successful in reproducing 



the observations. Dessart et al. found that [Siv] was relatively flat- 
topped in appearance, whereas the neon lines produce relatively 
symmetric and triple-peaked profile shapes. Ignace et al. found that 
the calcium line was quite asymmetric, with a strong blueshifted 
peak. It is possible to reproduce the gross profile morphologies for 
calcium and sulfur, but nothing like a triple-horned shape results 
for the neon lines at any phase. 

In order to match the profile shapes of the neon lines, an addi- 
tional dependence on radius or direction is required. Although our 
treatment of the bow shock head is very approximate, its shortcom- 
ings are not likely to suppress a three-horned morphology. Wolf- 
Rainer Hamann (private comm.) indicates that in the Potsdam mod- 
els of WR winds, the ionization balance is still not 'frozen-in' for 
WC models even out to the edge of their computational grid at 
200i?*. An example of this can be seen in Figure 7 of Grafener, 
Koesterke, & Hamann (2002) where at the lowest wind densities 
in the model, ion fractions are still changing with radius. We thus 
choose to introduce an additional scaling parameter associated with 
varying ionization as given by 

Q(r) oc (n c /n , c ) p , (34) 
for some constant exponent p. Given that n c oc r -2 , then Q oc 
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Figure 8. Similar to Fig.[7]except with application to the short period binary WR 11. The set of profiles at left assumes constant a constant ionization factor 
Q(r) = Qo; the set at right allows for non-constant Q(r) that best reflects the observed profile shapes (see text for details). The dashed line corresponds to 
an orbital phase of 150°, the model closest to 143° as determined for the time that the ISO data were obtained. 



Table 4. Wind and Binary Parameters for 7 Vel 



WR 11° = 7 Vel 


WC8+07.5 


M 


9 X 10~ 6 M yr" 1 


Voo 


1550 km s" 1 


ii% 


3.2 R e 


Me 


4.5 


Del 


10 


™0,e 


7.8 x 10 11 cm" 3 


e 


1/33 


Orbit 6 


d 


336 pc 


P 


78.53 days 


a 


1.2 AU 


e 


0.334 


i 


65.5° 


Interaction Region 




0.28 


/3 d 


85° 




750 i?*/rad 


r w /r c (CalV) 


4.2 


r-w/r-c (Nell) 


0.68 


r w /r c (Nem) 


0.36 


r w /r c (SIV) 


0.13 



a de Marco et al. 2000 b North et al. 2007 c T = 50120.4 MJD, and ISO 
data of 7 Vel was obtained on T = 50220.9836 MJD, hence the O star 
was at orbital phase 0.28. d Henley, Stevens, & Pittard 2005 



The right side of Figure[8]shows model profiles with p — 2 for 
Ca 3+ , 1 for Ne 2+ , 2 for Ne 1+ , and for S 3+ . Note that changing 



p to increasingly positive values tends to move the line formation 
deeper into the wind, because the line luminosity scales as 

L„oc J j v {r,$)Q{r)dV. (35) 

In addition, the larger density in the compressed layer, by a factor of 
4, tends to enhance recombination, and so for our models, we allow 
p — > p + 1 in the shocked layer. Comparisons to detailed radiative 
transfer calculations, such as those of the Potsdam group or CMF- 
GEN (Hillier & Miller 1999), are needed to assess the plausibility 
of the selected p values. Although motivated by the ionization con- 
siderations, Q(r) mainly acts as an additional weighting function 
to shift the predominant region of line formation in radius. 

Ultimately, our profile modeling still falls short of accurately 
reproducing the observed emission lines. The [CalV] line does be- 
come more extreme in its lopsided appearance for some phases. 
The [Nell] line can have a triple-horned appearance at some phases, 
which is a good sign; unfortunately, the triple-peaks are never as 
equal looking in appearance as observed. With Q a constant, the 
sulfur line continues to display mildly sloped or a mainly flat- 
topped appearance for restricted phases. 

There are two main problems with the modeling for 7 Vel. 
First, [Nelll] always shows a more or less double-horned morphol- 
ogy, never triple-horned as observed. Second, using the ephemeris 
of North et al. (2007), we have determined the orbital phase of 
7 Vel for the ISO line data. According to North et al., the O star 
orbits counter-clockwise on the sky, and so the CWIR spirals out- 
ward in a clockwise fashion. According to the ISO archive, the data 
were obtained on 1996 May 17 at 23:36:24, or MJD 50220.9836. 
Using the observation of North et al. on MJD 50120.4 and a period 
of P — 78.53 days, the O star is 28% of an orbit past periastron in 
time, which corresponds to an orbital phase of 143° . In our simula- 
tions this is close to a phase of 150° that is plotted as a dashed line. 
Thus, the second problem is that we cannot simulaneously get all 
four profiles with the overall grossly correct morphologies at quite 
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the same phase. Note that at a nearby phase of 120°, the [Nell] line 
does show a triple peak. 



4 CONCLUSION 

Our study concerns model calculations for forbidden emission line 
profile shapes that form in a colliding wind system. The approach 
recognizes that spherically symmetric winds produce optically thin 
and flat-topped emission line shapes produced in the constant ex- 
pansion flow. Deviations from that flat-top morphology represents 
an opportunity to infer information about the orbital properties of 
the binary and the colliding wind bow shock. 

Our approach makes a number of significant simplifications 
to make the problem tractable: the OB companion wind makes no 
contribution to the line emission, flow in the compressed layer is 
radial and with the same speed as the WR wind, the extent of the 
shocked layer is only approximate, the geometry of the bow shock 
head is treated poorly, the bow shock is assumed axisymmetric 
about the line of centers for the two stars, and more input regarding 
the ionization balance in the large scale wind is needed to model 
line shapes properly. 

In light of all of these shortcomings, the models do produce a 
remarkable degree of diverse profile morphologies. Although even 
qualitative matches to observed forbidden lines from WR 147 and 
7 Vel are far from satisfactory, the attempt does appear to capture 
a number of trends and shows potential as a tool for deducing or 
limiting orbital and wind parameters. 

Although of lower spectral resolution than ISO, we deter- 
mined that the Spitzer IRS could probably detect modulations of 
IR forbidden line shapes with orbital phase. Unfortunately, the IRS 
will not be available during Spitzer 's "warm" cycle. It is expected 
that the James Webb Space Telescope will have a mid-infrared 
spectrograprQ with a resolving power similar to ISO. In addition, 
the [Nell] emission line is observable from ground-based observa- 
tories, as for example in the study of Smith & Houck (2001). In par- 
ticular, their sample included WR 146 which has quite broad lines 
owing to a wind terminal speed of nearly 3000 km s" 1 . So there is 
the possibility that future observations will provide new high qual- 
ity data of these and other systems for which our diagnostics will 
be relevant. 
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